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PREFACE 


The  techniques  for  computing  confidence  intervals  for  spectral  estimates 
from  the  Coastal  Data  Information  Program  were  developed  at  the  University  of 
Wyoming  during  the  period  from  August  1982  to  August  1983.  This  report  was 
done  at  the  request  of  and  was  funded  by  the  US  Army  Engineer  Division,  South 
Pacific,  as  part  of  the  Coastal  Field  Data  Collection  Program  at  the  Coastal 
Engineering  Research  Center  (CERC)  of  the  US  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES).  The  work  was  performed  by  Drs.  Michael  E.  Andrew,  Re¬ 
search  Statistician,  CERC,  and  Leon  E.  Borgman,  Professor  of  Statistics, 
University  of  Wyoming. 

This  report  was  revised  by  Dr.  Andrew  and  edited  for  publication  at 
CERC  under  the  general  supervision  of  Drs.  Dennis  R.  Smith,  Chief,  Prototype 
Measurement  and  Analysis  Branch,  William  L.  Wood,  Chief,  Engineering  Devel¬ 
opment  Division,  and  Dr.  Robert  W.  Whalin,  Chief,  CERC. 

Commanders  and  Directors  of  WES  during  the  conduct  of  this  study  and 
preparation  and  publication  of  this  report  were  COL  Tilford  C.  Creel,  CE,  and 
COL  Robert  C.  Lee,  CE.  Technical  Director  was  Mr.  Frederick  R.  Brown. 
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METHODS  FOR  COMPUTING  CONFIDENCE  INTERVALS  FOR 
SPECTRAL  ESTIMATES  IN  MONTHLY  REPORTS  OF  THE 
CALIFORNIA  COASTAL  DATA  INFORMATION  PROGRAM 

PART  1:  INTRODUCTION 

T’n 2  limits  containing  a  certain  quantity  with  a  probability  of  1-6 
are  called  the  (l-«)  •  100  percent  confidence  limits  for  that  specific 
quantity.  The  interval  between  the  confidence  limits  is  called  a  confidence 
interval.  The  purpose  of  this  paper  is  to  present  methods  for  finding  con¬ 
fidence  intervals  for  the  quantities  estimated  in  the  monthly  reports  of 
the  Coastal  Data  Information  Program  sponsored  by  the  US  Army  Corps  of 
Engineers  and  the  State  of  California  Department  of  Boating  and  Waterways.  - 
The  theory  presented  here  pertaining  to  confidence  limits  for  estimates  of 
the  spectral  energy  has  been  in  existence  for  some  time  and  has  been  veri¬ 
fied.  A  good  discussion  of  its  validity  appears  in  a  recent  paper  by 
Donelan  and  Pierson  (1983).  The  theory  concerning  the  probability  law  for 
the  longshore  component  of  radiation  stress  estimates  and  associated  con¬ 
fidence  intervals  (as  well  as  the  approximate  confidence  intervals)  is 
relatively  new  and  has  not  been  widely  discussed  in  any  applied  sense. 
However,  some  theoretical  results  have  been  presented  by  N.  R.  Goodman 
(1957)  on  the  probability  law  for  estimates  of  the  cospectrum. 
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PART  2:  CONFIDENCE  INTERVAL  FOR 
PERCENT  ENERGY  ESTIMATES 


It  has  been  shown  that,  under  the  assumption  of  a  stationary  Gaussian 
sea  surface,  the  normalized  spectral  density  estimates  obtained  by  finite 
Fourier  techniques  and  frequency  smoothing  have  a  chi-square  divided  by 
degrees  of  freedom  probability  distribution  (Borgman  1973).  Furthermore, 
this  theory  has  recently  been  demonstrated  to  yield  accurate  estimates  of 
the  sampling  variability  for  the  spectral  estimates  (Donelan  and  Pierson 
1983). 

If  S(f)  denotes  the  spectral  estimate  from  an  average  of  k/2  raw  spectral 
lines  centered  at  frequency  f,  then 

Iffy  *  t  (2'1) 

where  x£  represents  a  chi-square  random  variable  with  k  degrees  of 
freedom.  Thus 


Pr^xk ,a/2  <  SlfT  <  Xk ,1 -a/2^  =  1_a 


where  y}  is  the  value  of  the  chi-square  variable  for  which 

'  K  y  9 


Prt4  <  <,a]  =  a 


(2.2) 


(2.3) 


The  estimates  provided  in  the  reports  are  in  terms  of  percent  of  the 

A 

total  energy.  Thus  it  is  necessary  to  find  a  relation  between  S(f)  and 
PE(f),  the  percent  of  total  energy  for  the  frequency  interval  centered  at  f. 


PE(f )  =  (Af  j  Sj  M 
m 


(2.4) 
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where 


Af  -  digital  increment  in  the  frequency  domain 
=  raw  spectral  line  for  frequency  mAf 
TE  =  total  energy  in  spectrum 

Note  that  the  summation  over  m  refers  only  to  that  frequency  interval 
centered  at  f.  Thus, 


(2.5) 


PE(f)  =  |  Af  S(f )  -14j| 


(2.6) 


since 


s<f>  -  t  i  sm 

m 


PrrAfl00Y2  .WlAfM  ]  = 

KrL  2  TE  xk,a/2  ^[fJ  2  TE  xk,l-ct/2J 


(2.7) 


(2.8) 


k 


u 


rPE(f]_  2TE  <  s(f,  <  PE(f)  2TE  =  , 

C4fl005[^  S(f)  "IOOxJ^J 


(2.9) 


PE(f)  2TE  PE(f )  2TE 

Af  100  CC/2  ’  10°  xk ,a/2 


(2.10) 


yields  a  (1-ct)  x  100  percent  confidence  interval  for  S(f).  However,  since 
the  reports  are  in  terms  of  percent  of  total  energy,  consider  the  following 


expression: 


,  k  ;  kAf  S(f) 

Cl ^2  2  TE 


Pr[PE( f )  -r-5 -  <  ^  x  x  100  <  PE(f) 


xk,a/2 


]  =  (2.11) 
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Then 


PE(f)  ~7 - 

Xk,l -a/2 


k  ,  PE(f)  t~t~  k 


k  ,a/2 


(2.12) 


is  a  (1-a)  •  100  percent  confidence  interval  for  the  quantity 


MiLii  ioo 

2  TE  1UU 


(2.13) 


The  quanitity 


PE (f )  =  ~  •  100 


(2.14) 


represents  the  discrete  approximation  of  the  energy  in  the  frequency  interval 
centered  at  f  containing  k/2  raw  spectral  lines. 

The  monthly  reports  for  the  Coastal  Data  Information  Program  contain 
nine  spectral  estimates  for  each  location  and  time.  The  spectral  estimates 
are  given  by  period  bands  and  result  from  varying  values  of  k,  as  shown 
in  Table  1.  Also  listed  are  the  values  for  and  xj<  ]-a/2  for  the 

appropriate  values  of  k  by  period  band. 


Table  1 

Coastal  Data  Information  Program  Spectral  Estimates 


Period 

k 

xk,.05 

xk,  .95 

22+ 

92 

71.76023 

116.50839 

22-18 

20 

10.84948 

31.41630 

18-16 

16 

7.96032 

26.30112 

16-14 

18 

9.38799 

28.87667 

14-12 

24 

13.84490 

36.42149 

12-10 

34 

21.66101 

48.61088 

10-  8 

52 

36.43630 

69.82803 

8-  6 

84 

63.87607 

106.39197 

6-  4 

172 

142.6713 

203.59982 
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Table  2  lists  the  estimates  from  the  Crescent  City  array  for  May  1,  1981 
at  time  2120.  Also  listed  are  the  lower  and  upper  90%  confidence  bounds  for 
the  quantity  PE  .(see  Appendix  B  for  a  copy  of  the  original  pages  from  the 
report) . 


Table  2 

Crescent  City  Spectral  Estimates 


Period 

A 

PE 

Lower 

Upper 

22+ 

3.0 

2.369 

3.846 

22-18 

0.6 

.382 

1.106 

18-16 

0.6 

.365 

1.206 

16-14 

3.0 

1.870 

5.752 

14-12 

43.8 

28.862 

75.927 

12-10 

20.7 

14.478 

32.492 

10-  8 

12.3 

9.160 

17.554 

8-  6 

7.4 

5.842 

9.731 

6-  4 

9.1 

7.688 

10.971 

Using  the  values  in  Table  1  and  Equation  (1.2)  to  obtain  confidence 


intervals  on  the  quantity  defined  as  PE  is  relatively  simple  for  any  set  of 
estimates  taken  from  the  monthly  reports. 


PART  3:  PROBABILITY  LAW  FOR  THE  LONGSHORE 
COMPONENT  OF  RADIATION  STRESS 


The  longshore  component  of  radiation  stress  denoted  by  S^y  has  been  shown 
to  be  a  function  of  the  cospectrum  of  the  offshore  and  longshore  surface 
slopes  (Seymour  and  Higgins  1978).  Assuming  linear  wave  theory,  it  is 
possible  to  determine  the  probability  law  for  Sxy  in  a  closed  mathematical 
form.  Denote  the  two  surface  slope  time  series  by 


n  (nAt)  =  offshore  slope 
n  (nAt)  =  longshore  slope 


(3.1) 


U  (mAf)  -  iV  (mAf)  =  At  [  n  (nAt)  exp(-i2Trmn/N) 
x  x  n=0  x 


(3.2) 


U  (mAf)  -  iV  (mAf)  =  At  l  n  (nAt)  exp( -i 2mmn/N) 

y  y  n=o  y 

which  are  the  discrete  Fourier  transforms  of  the  two  slope  components, 
estimate  for  the  cospectrum  of  with  r;y  is 


Cxy(mAf)  =  -A- 


Uv(mAf)  U  (mAf )  +  Vv(mAf)  V>;  (mAf)) 


(3.3) 


An  estimate  for  S  as  a  function  of  frequency  is 
xy 

Sxy(mAf)  =  P3-k^  Cxy(mAf) 


where 


k  =  wave  number 

p  =  water  density 

g  =  acceleration  due  to  gravity 

n(k)  ‘  2  Sinh  (2k h) 
h  =  measured  average  depth 


(3.5) 
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and  n(k)  and  k  are  functions  of  frequency  according  to  the  followin';  relation 


(2  :irr  f )  =  gk  tanh(kh) 


.( 


It  can  be  shown  that  Ux(nrf),  U  (irrf),  Vx(nrf),  V^(in  f)  are  distributed 
multivariate  gaussian  with  mean  zero  and  variance  covariance  matrix: 


S  (m  f) 

XX  v  ' 

C  (nrf) 
xy v  ‘ 

0 

0 

C  ( nr  f ) 
xy v  ' 

S  (in  f) 

yy 

0 

0 

0 

0 

Sxx(m  f) 

C  (m  f) 
xyv 

0 

0 

S  (m  f) 

yy 

where 

S(f)  =  The  spectral  density  of  the  sea  surface 
5  (f)  =  The  spectral  density  of  the  offshore  slope 

X  X 

S  (f)  =  The  spectral  density  of  the  longshore  slope 

Also  the  four  random  variables  are  independent  for  different  values  of  m  in 
the  interval  0  •-  m  <  N/2. 

In  order  to  obtain  the  probability  law  for  C  (nrf)  it  is  necessary  to 
form  the  characteristic  function.  Define  the  characteristic  function  for  a 
random  variable  zl  to  be 


A2,(t)  '  i[eitzl] 


stands  for  statistical  expectation  and  i  •/ - 1.0  . 

vim  f 


'hi is  if  wo  lot 


n  <■  1 1 


1 1 


I  L. 
\ 


FT' 


i  tzl 


r  /exp  ji  t[u  ii’m  f  ;  I !  ( in '  t )  ♦  V  (nrf)  V.Oirf)]^ 


Then  let 


cos-  sin- 
cos-  sin-1  ’ 


(4.22) 


and  then 


Sxy(f)  E(f)  _  S(f) 

Sxpl  ETfT  '  "  Slf) 


(4.23) 


It  has  already  been  stated  that  S(f)/S(f)  has  a  chi-square  divided  by  degrees 

of  freedom  probability  distribution.  Then  for  S  (f)  >  0 

xy 


Pr[ 


k,  -72 
k 


s  f 
s  14 

xyv 


’k  ,1-  72  q  ..  , 

“  _  — j  1  -  t 


(4.24) 


or 


k  S  (f) 
Pr[  ---V.-jL. 

’k ,1-  72 


k  S  (f) 

S  ,(f)  <  ]  •  1-a 

r  xk,r72 


xy 


(4.25) 


Note  that  if  S  (f)  •  0.0  then  the  inequalities  are  reversed.  Thus 
the  resulting  confidence  bounds  are  reversed. 

The  expression  for  can  be  rewritten  by  letting  A  =  r-  -  f>  and 


sin  2- 
s  i  n  2 

sin  2 

sin  2 (  + ‘ ) 

sin  _ ... _  (4.26) 

sin  2  cos  2'  +  cos  2-  sin  2\ 

“  1/  * 


wnnm  -  cos  3  *■  sin  ?'  cot(?  ).  Then  the  confidence  interval  can  be 

written  as  follows: 


J 

•J 


•J 
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several  different  values  of  the  frequency  averaging  bandwidth  then  plotting 
the  coherence  estimate  vs.  bandwidth.  The  coherence  estimate  will  be  equal  to 
one  for  one  spectral  line  per  averaging  band.  It  will  then  drop  off  rapidly 
toward  zero  for  increasing  bandwidth  until  it  reaches  a  plateau.  Beyond  this 
plateau  it  will  again  drop  off  to  zero  with  increasing  number  of  lines  per 
average.  This  is  due  to  the  smoothing  across  features  in  the  true  spectral 
densities.  The  largest  bandwidth  that  will  keep  the  coherence  on  the  plateau 
is  a  good  choice  since  coh(f)  is  stabilizing  around  its  theoretical  value. 

Also, the  largest  bandwidth  on  the  plateau  provides  the  highest  degrees  of 
freedom  for  the  estimates  (Borgman  1973).  Estimates  for  coh(f)  are  not  avail¬ 
able  in  the  monthly  reports.  Due  to  this  difficulty  it  is  not  possible  to 
obtain  approximate  values  for  t-j  and  t^.  However  it  is  possible,  under  certain 

assumptions,  to  obtain  rough  confidence  intervals  for  S  (f).  The  monthly 

xy 

reports  contain  values  for  the  apparent  angle  n  which  is  related  to  S  (f) 

xy 

by  the  following  equation: 

S  (f)  =  E(f)  n(f )  cosO  sinO  (4.19) 

where  E(f)  is  the  energy  contained  in  the  frequency  interval  centered  at  f 
(Seymour  and  Higgins  1977). 

Let  E[f]  be  approximated  by 

E(f)  =  S(f)  6(f)  (4.20) 

where  -(f)  is  the  width  of  the  frequency  interval  centered  at  f  ,  in  this  case 
6(f)  =af  k/2.  Also  let  E[f]  be  estimated  by 

E(f)  =  S(f)  6(f)  (4.21) 


> 
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However,  the  function  f*(z)  contains  an  unknown  parameter  D  where 


D-  = 


S  (f)  S  (f)  -  C2  (f) 

xx v  '  yy  xyv 

C2  (f) 
xy  ' 

S  (f)  S  (f) 

xx  _  1-0 

C2  (f) 
xy  ' 


-  coh(f) 


- 1 


1.0 


(4.15) 


The  expression  coh(f)  stands  for  the  coherence  of  the  x  slope  with  the  y  slope. 
Coherence  is  defined  to  be 


coh(f ) 


C2  (f)  +  q2  (f) 

S  (f)'  S  rh 
xx  ‘  yy 


or  the  modulus  squared  of  the  cross  spectrum  divided  by 
densities.  In  the  case  of  x  slope  and  y  slope  the  quad 
is  theoretically  zero,  so 


(4.16) 


the  two  spectral 
spectrum  q  (f) 


coh( f ) 


C2  (f) 
xy 


S  Jf) 

XX 


s  “TfT 
yy 


(4.17) 

2 


Since  the  coherence  is  not  known,  it  is  not  possible  to  determine  D  and 
thus  t-j  and  t2  exactly.  If  t-j  and  t2  are  found  to  satisfy  the  two  conditions 
using  an  estimated  value  of  coh(f),  it  is  not  certain  that  the  confidence 
interval  will  then  have  the  specified  confidence  coefficient  (1-0  •  100 
percent.  The  estimated  coherence  has  the  following  form: 


coM  f '• 


S  (f) 
x  .< 


S  (f) 
yy 


(4.18) 


To  "'ini  ice  an/  sources  of  error  due  to  the  use  of  an  estimated  coherence  in 
obtaining  /a Urns  for  t. ^  and  ,  it  is  desirable  to  have  the  best  possible 
estimate  for  coh(f).  This  is  done  by  computing  the  coherence  estimate  for 
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Case  2:  Assume  <  0  ,  then  the  result  is  the  same  as  in  Case  1. 

A 

Case  3:  Assume  t,  < 0  ,  t9  >  0  >  and  S  (f)  >  0  ,  then 
i  xy 

^  (f) 

Pr[tl  1  V?!"  T>]  ' 

s„v(f>  i  „(f) 

*  Pr[t'  i  s^ITT  ■=  0]  +  Pr[°is^r?Ti  V 

s  (f)  S  (f) 

=  PrC-  <  S  (f)  <  -Sf - ]  +  Pr[^ -  <  S¥  (f)  <  co]  (4.10) 

Ay  t-|  1*2  xy 

=  l-x 

Then  for  t-j  0,  t?  ^  0  and  S  (f)  _>  0,  (1-a)  •  100  percent  confidence 
interval  is 

S  (f)  S  (f) 

(-■,  )  and  - ,  «)  (4.11) 

Z]  Z2 

Case  4:  If  t-|  0,  t^  ■  0  and  Sxy(f)  <  0  ,  then  the  intervals  become 

S  (f)  S  (f) 

(-,  -X{-  )  and  (-x{— ,  -)  (4-12) 

Z2  n 

The  na’or  difficulty  in  using  the  above  intervals  is  in  the  determination 


of  the  values  for  t. .  and  t~.  The  conditions  for  this  are 


_  r  w  . 

f*(z)  '=  27  !  Qz(t)  e'ltzdz 


_  -  f  (A) 

—I  \  ® 


00  -itz 


00  {t-ia}v(t-ib)v 


2rr  'O'  j 


/ v \ 2v  -  /  \ 

=  (p)  f(z) 


0  az  v-1  ,  .  , v ,  v-m-1  , 

_  /V\2v  e  y  (v+m-1)!  z  1 

'rv  7771  1  1  L_  7u-m-l  1 1  ml 


D  Tv^irr  m^0  Cv-nT-1)!  m!  (a+b)v+m 


(4.6) 


a  =  a  ,  z  <  0 

a  =  b  ,  z  >  0 


having  obtained  f*(z)  it  is  possible  to  find  t^  and  such  that 

S  (f) 

Pr[tl  -  S^TfJ  -  f'2]  =  1-ct 


(4.7) 


In  order  to  invert  inequality,  it  is  necessary  to  break  the  values  of 
t-j  and  t2  into  cases. 


Case  1  :  Assume  t-j  0  ,  then 

7(0  s  (f> 


(4.8) 


Then  the  (1-0  •  100  percent  confidence  interval  for  S  (f)  is 

*y 

IS  (f)  S  (f)l 

j  _xy__.  *y_ .  . , 

*-  *  f 

1  7  7  j 


(4.9) 
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PART  4:  A  CONFIDENCE  INTERVAL  FOR  THE  LONGSHORE 
COMPONENT  OF  RADIATION  STRESS 


In  the  section  of  this  report  where  the  probability  law  for  the  longshore 

component  of  radiation  stress  was  developed,  the  random  variable  z,  where 

z  =  v  S  (f)/3(f),  was  obtained.  Here  consider  the  variable 
xy 

2  ■  Vf>'Vf>  i4-1’ 

where  S  (f)  is  the  estimate  defined  in  the  previous  section,  and  Sx^(f)  =  »(f 
is  the  theoretical  value.  The  characteristic  function  for  the  new  variable 
z  is  from  equation  2.28. 


Qz(t)  =  (1 


it  ,  t2D21-v 
+  — 2—  ) 


(4.2) 


where 


D?  = 


s..(0  svv(f)  -  c;y(f) 


The  function  Q  (t)  can  be  rewritten 


Qz(t)  =  (^-)2v  { t-i a>_v{ t-i b> 


-v 


where 


(4.3) 


(4.4) 


3  =  25  ([)  +  A/D*  +  *) 

,  (4.5) 

b  =  ^  (l  -  rW  +  4) 

Since  (f)  v  is  constant,  it  is  convenient  to  obtain  the  probability  density  of 
z  from  the  results  in  the  previous  section.  Thus 


If  for  the  curve  r  we  choose  t  =  Re  0  £  o  ^  jt,  then  the  conditions  for 
using  the  residue  theorem  are  fulfilled  (see  Appendix  A  for  proof),  and 


e~itz  q  (t)dt  =  -2-rri  a_-j 


(3.4^ 


where 


a-l  =  1ini  7^TTT  ~T\  Ut-ib)V1tz  Q  (t)> 
t+ib  v  at 

Since  Q  (t)  =  (t-ia)~v(t-ib)*v 


(3.41 


a-l  =  lim  f^-Tyi  ~Vr  {e"1tZ  (t-ia)~v)  (3-4' 

t->ib  v  ’’  at  ■ 

The  derivative  can  be  evaluated  using  the  formula  in  equation  (3.47). 

For  any  two  functions  f(x)  and  g(x)  such  that  the  required  higher  order 
derivatives  exist,  it  is  possible  to  show  that  (Abramowitz  and  Stegun  1970) 
~n  f  „  n  „  ~,n-k  f  -  k 

9  f'9  _  v  /n^  9 _ I  9g  (3  4 

n  n-k  k 

3xn  k=0  3xn  K  3xK 


3V  ^  e  ^z  (t-ib)  ^  ^v-1 ) ( _-j z^v-m-l  e~ i tz 

3tV_1  m=0  m 


,  , sm  (v+m-1 ) !  ,  ^  -u 
('1}  “TvTjT  (t_  b 


?bz  V  (v+m- 1 ) !  (-l)v+m  4  zv'f"~1 


[v-m-1 ) !  m! 


(b-a) 


+(v+m) 


Then  for  z  0 


f(t)  =  -  2-1  a., 


ebz  (v+m-1 ) !  z'-’-"’-' 


e Q  (t)dt  =  I  e  Q(t)dt  +  e',tz  Q(t)dt  (3.38 

c  z  Jr  z  >1  z 


where  F  is  the ‘semi circle  part  of  c  and  L  is  the  line  segment.  Now  if  it 


can  be  shown  that 


lim  e'1tz  Q  (t)dt  =  0 

R->oo  Jf 


lim  e"1tzQ  (t)dt  = 


(3.39 


R->co  >  Q 


Qz(t)dt 


(3.40 


However,  by  the  Residue  Theorem  of  complex  analysis  (Mardsen  1973),  the  inte 


(3.41 


lim  e'1  1  Q  (t)dt  =  +  2tt1  a_- 
R-**  >  c 


where  a.]  is  the  residue  of  the  function  e"1tz  Qz(t)  at  the  singularity 
contained  in  the  curve  c.  The  function  e~1tz  Qz(t)  has  a  singularity  at 
ia  and  one  at  ib.  Both  of  these  points  are  said  to  be  poles  of  order  v  . 
Thus  the  function  e”ltz  Qz(t)  has  residues  equal  to  (Spiegel  1964) 


Jit  T^T)!  3tv_1  6  UZ  Qz(t)} 


(3.42 


where  t^  =  ia  and  ib.  If  the  residue  is  in  the  upper  half  of  the  plane,  then 
the  integral  has  value  +2  i  a  j.  If  it  is  in  the  lower  half,  then  the  inte¬ 
gral  is  - 2  ’ i  a  i 

It  is  necessary  to  consider  two  cases  in  order  to  compute  the  integral 


If  the  formula  is  factored  inside  of  the  brackets,  the  function  becomes 


Q,(t)  =  (t-ia)'v(t-ib)~v 


(3.32) 


where 


(3.33) 


(3.34) 


According  to  the  inversion  theorem  for  characterise c  function  (Rao,  1973),  the 


probability  density  function  of  the  variable  z  =  vS  (f)/S(f)  is 

xy 

/*  oo 

f(z)  =  27  e'UZ  Qz(t)dt 


=  07  e”itz  (t-ia)~v(t-ib)~v  dt 

-oo 


(3.35) 


(3.36) 


It  is  necessary  to  make  use  of  the  theory  of  contour  integrals  in  order  to 


solve  for  f(z).  Consider  the  integral  on  the  curve  c 


e'itz  Qz(t)dt 


(3.37) 


where  the  curve  c  is  either  of  the  semicircles  of  radius  R  centered  at  the 
origin  and  the  line  segment  from  -R  to  R,  as  indicated  below. 


bince  the  curve  c  is  closed  and  since  Q  (t)  is  analytic  everywhere  but  at  the 
points  ia  and  ib,  then 


B(mAf)  =  pg 


(3.24) 


(S  (mAf)  S  (mAf)  -  C  2(mAf)P 
k  '  xx  yy  xy  '  ' 


I 


Since  n-^-  C  (mAf)  is  independent  for  different  values  of  m  it  is  possible 
k  xy 

to  use  another  property  of  characteristic  function.  That  is,  if 


z  =  1  2m  (3.25! 

m 

where  z  is  independent  for  different  m,  then 

«2(t)  =  {Qz(t)}m  (3.26! 

Under  the  assumption  that  the  parameters  a(mAf)  and  [(mAf)  can  be  taken  to  be 
constant  over  the  interval  of  frequency  smoothing  if 


23  =  ®  V  c  (m4f)  , 
v  '  k2  xy'  '  ’ 
m 

then 

oz3{t)  =  n  -  itf-.  e  ^P-rv 


(3.27; 

(3.28) 


In  the  last  equation  <(f)  and  3(f)  denote  the  values  of  the  two  parameters  at 
the  midpoint  of  the  frequency  interval  of  smoothing. 

In  order  to  simplify  the  function  |>  ,  if  the  new  random  variable  z 
is  to  be  defined  as 


z  =  .S  (  f )/.-’(  f) 
xy 

then 


(3.29) 


?2(t)  ’  '1  •  u  4vj  ♦  t?A‘  (3. 

let.  .(f)-  (f)/  (f)  and  (f)  /  0.  Then  ,  is  a  function  of  only  two  param- 

L. 

eters  and  .(f)  that  is 


Q  (t)  -  ‘  1  -  i,(f)t  +  t:i- 


(3.31 


. . .  ■■  »i  »'u » 'j  ■[■  i.1  i»»"  »n 


m 


«■ . 


« 


Thus 


Qzl(t)  =  |I  -  2it  CM f 


-k 


=  {|C*|2}' 

=  |c*H 


(3.18) 


=  {1  -  it  Cxy(mAf)  +  ^  (Sxx(mAf)  Syy(mAf)  -  C^^mAf))}"1 

If  Sxy(mAf)  is  smoothed,  then  the  smoothed  estimate  will  have  the  following 
form: 


S  (f)  =  -  Y  S  (mAf) 
xyv  '  v  L  xy  ' 
m 


(3.19) 


where  the  sum  is  over  a  specified  frequency  band  with  v  spectral  lines.  Since 


S  (mAf)  is  a  function  of  C  (mAf),  then 
xy  xy 

s  (f)  =  fiR  V  nM  c  (n^f) 
xyv  '  v  £  k2  xy'  ' 

One  of  the  properties  of  characteristic  functions  is  that 


<W*>  *  «z<T> 


Thus  if  z2  =  pg  C  (mAf) ,  then 
k  xy 

n(k) 


(3.20) 


(3.21) 


and 


Qz2(t)  =  (1  -  it  pg  ^  Cxy(mAf) 

+  ^  P9  (Sxx^mAf)  Syy(tTlAf) 

=  (1  -  it  ot(mAf)  +  t282(mAf)}  1 


a(mAf)  =  pg  Cxy(mAf) 


Cxy2 (mAf) ) } _ 1 


(3.22) 


(3.23) 
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Cxy(mAf) 

S  (mAf) 

XX 

0 

0 

S  (•  *f) 
y.y 

C  (mAf ) 
xy 

0 

0 

0 

0 

C  (mAf) 
xy  ’ 

Sxx(“f) 

0 

0 

S  (mAf) 

yy 

Cxy(mAf) 

'  c*  o' 

I  -  2 i t  CM  =  ( 

,  c  c* 

1  -  i|  Cxy(m.lf)  -if  Sxx(m4f) 

C*  =  .  .  ( 

-if  Syy(«f)  I  -  if  cxy(«f) 

f  C*  0 

I I  -  2i t  CM  I  = 

0  C*j 

=  [C*|  |C*| 

=  |C*|2  (: 

|c*l  =  0  -  T  Cxy(mAf))-  +  SXX(mAf)  Syy(mAf) 

=  1  -  it  Cxy(mAf)  -  Cxy2(m;,f)  +  Sjm'.f)  Syy(mAf) 

=  1  -  it  C  ( m f )  +  (S  (ni'.f )  Sw  (mAf)  -  C  ?(m\f)) 


f 


where 


J 1'  =  [U  (mAf) ,  l)  (mAf ) ,  V  (mAf) ,  V  (mAf)] 

x  y  x  y 


i  -yrc"lu 


Since  f  (U)  =  |2TyT| C | h  e 


(3.10) 


Qzl(t) 


it  U"MU  1  ,  -%U*C_1U  ... 

! - l2^F|C|ls  e  dU 


1  ,  _-%(U'C-1U  -21t\TMU)  . 

TuV\t\2  - 


1  .  -4  U'[C_1-2it  M]U  du 

x^rici^6  - 


(3.11) 


Qzl(t) 


\-h  (  00  t  e 


-^r[C“l-2it  M]U 


( 2tt ) ^  |C_1+2it  M | 


‘+2it  Ml  2 


ICl-*4  [  C~  1  -  2  i  t  m;"'2 


I  -  2it  CM | " 5 


(3.12) 
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r*k  S  (f)  *  k  S  (f)l 

xy '  xy  ' 

xk,l-,/2  vk,,/2  j 


(4.27) 


It  is  possible  now  to  compute  confidence  intervals  for  S  (f)  given  arbitrary 

xy 

choice  of  A.  If  A  is  assumed  to  be  small  in  absolute  value  compared  to  cot(2e), 
then  -*  will  be  close  to  1.0.  However,  for  u  close  to  zero  cot(2f)  becomes 
very  large  and  can  dominate  the  value  of  •  *.  In  this  case  it  may  be  desirable 
to  choose  a  value  for  *  that  is  different  from  1.0,  according  to  an  arbitrary 
choice  for  A.  Table  3  lists  some  values  of  *  for  typical  values.  For 
negative  values  simply  look  under 


Table  3 

Values  of  l*  for  Typical  o  Values 


3  ' 

7 

9 

11 

_ 15_ _ 

47 

-2.0 

.33 

.72 

.78 

.82 

.88 

1.00 

-1.0 

.67 

.86 

.89 

.91 

.94 

1.00 

-  .5 

.83 

.93 

.95 

.96 

.97 

.99 

0 

1.0 

1.0 

1 .0 

1  .0 

1.0 

1.0 

.5 

1.17 

1.07 

1.05 

1  .04 

1.03 

1.01 

1.0 

1.33 

1.14 

1.11 

1.09 

1.06 

1  .00 

2.0 

1.66 

1.28 

1.21 

1.17 

1.12 

.99 

The  following  estimates  were  taken  from  the  Crescent  City  array  on  May  1, 
1981  (time  2120)  (Appendix  B).  Table  4  lists  the  period  bands,  associated 
midpoint  frequencies,  wave  numbers  and  value  for  the  function  n(f).  Table  5 
lists  the  reported  angles  measured  clockwise  from  true  north,  the  associated 
m  values  that  are  counter-clockwise  from  beach  normal  (positive  x  axis),  and 
the  S  values  computed  with  the  following  formula: 

S  (f)  =  [ff)  n(f)  cos'1  sin"  (4.28) 

xy 
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Table  4 

Midpoint  Frequencies  and  Wave  Number  Values 


Per  i  od 

Midpoint  Frequency 

Wave  Number 

n(f). 

22  + 

.00273 

.0020:' 

.51580 

22-13 

.05050 

.01027 

.99797 

13-16 

.05903 

.01 104 

.99623 

16-14 

.06696 

.01306 

.99377 

14-12 

.07738 

.02412 

.98897 

12-10 

.09167 

.03385 

.97860 

10-  8 

.11250 

.05098 

.95325 

8-  6 

. 1 4580 

.08567 

.88238 

6-  4 

.20833 

.17484 

.68726 

Table  5 

Reported  and  Adjusted  Directional  Estimates 


Peri od 

Reported  Angle* 

0 

S  (f) 
x.y '  ‘ 

22  + 

22-18 

50.9 

-9.1 

-1  .048 

18-16 

15.0 

-45.0 

-3.350 

16-14 

44.4 

-15.6 

-8.655 

14-  12 

66.9 

6.9 

57.903 

12-10 

32. 0 

12.0 

46.173 

10-  8 

65.9 

5.9 

13.437 

8-  6 

64.5 

4.5 

5.724 

6-  4 

57.4 

-2.6 

-3.176 

*  Angles  in  degrees. 

Table  6  lists  the  b  (f)  values  from  Table  5  together  with  the  lower  and 
upper  confidence  bounds  for  varying  guesses  of  '  .  The  chi-square  values 
and  degrees  of  freedom  are  the  same  as  for  the  spectral  estimates  (Table  1). 
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Table  6 

Confidence  Intervals  for  S 


Peri  od 

S 

xy 

★ 

Lower 

Upper 

22-18 

-  1.048 

1.22 

-  2.36 

-  .82 

18-16 

-  3.350 

1.03 

-  6.93 

-  2.10 

16-14 

-  8.655 

1.12 

-18.58 

-  6.04 

14-12 

57.903 

.71 

27.09 

71.24 

12-10 

46.173 

.83 

26.80 

60.16 

10-  8 

13.437 

.66 

6.60 

12.66 

8-  6 

5.724 

.56 

2.53 

4.21 

6-  4 

-  3.176 

1.77 

-  6.78 

-  4.75 

22-18 

-  1.048 

1.11 

-  2.14 

-  .74 

18-16 

-  3.350 

1.02 

-  6.87 

-  2.08 

16-14 

-  8.655 

1.06 

-17.59 

-  5.72 

14-12 

57.903 

.86 

32.15 

86.29 

12-10 

46.173 

.92 

29.71 

66.68 

10-  8 

13.437 

.83 

8.31 

15.91 

8-  6 

5.724 

.77 

3.48 

5.80 

6-  4 

-  3.176 

1.38 

-  5.28 

-  3.70 

22-18 

-  1.048 

1.0 

-  1.93 

-  .67 

18-16 

-  3.350 

1.0 

-  6.73 

-  2.04 

16-14 

-  8.655 

1.0 

-16.59 

-  5.39 

14-12 

57.903 

1.0 

38.16 

100.34 

12-10 

46.173 

1.0 

32.29 

72.48 

10-  8 

13.437 

1.0 

10.01 

19.17 

8-  6 

5.724 

1.0 

4.52 

7.53 

6-  4 

-  3.176 

1.0 

-  3.83 

-  2.68 

22-18 

-  1.048 

.89 

-  1.72 

-  .60 

18-16 

-  3.350 

.98 

-  6.59 

-  2.00 

16-14 

-  8.655 

.94 

-15.59 

-  5.07 

14-12 

57.903 

1.14 

43.50 

114.39 

12-10 

46.173 

1.08 

34.87 

78.28 

10-  8 

13.437 

1.17 

11.71 

22.43 

8-  6 

5.724 

1.22 

5.51 

9.19 

6-  4 

-  3.176 

.62 

-  2.37 

-  1.66 

22-18 

-  1.048 

.78 

-  1.50 

-  .52 

18-16 

-  3.350 

.96 

-  6.46 

-  1.96 

16-14 

-  8.655 

.87 

-14.43 

-  4.69 

14-12 

57.903 

1.29 

49.22 

129.44 

12-10 

46.173 

1.16 

37.46 

84.08 

10-  8 

13.437 

1.34 

13.41 

25.69 

As  apparent,  the  intervals  can  be  very  unstable  for  changing  values  of 

\  .  Some  of  the  intervals  do  not  even  contain  S  (f).  This  demonstrates 

xy 

that  great  care  must  be  taken  in  the  i nterpretation  of  such  intervals. 

However,  if  one  is  willing  to  assume  a  value  for  the  apparent  angle  ■■ 
and  thus  a  value  for  \  ,  then  this  method  is  straightforward  and  relies  only 
on  the  well-known  chi-square  random  variable. 

If  one  is  willing  to  assume  a  value  for  the  coherence  that  is  consistent 
with  a  specific  model  for  the  directional  spectrum,  then  it  is  possible  to 
compute  the  constants  and  t  and  obtain  the  confidence  interval.  The 
next  section  of  this  report  will  present  such  an  approach. 


1'ART  S:  CONK  I  DI-NCI-!  !  NTl-RVAl.S  FOR  Till-.  l.ONOSItORF 
COMPONKNT  OF  RADIATION  SI  RTSS  IIV  Till. 
U'RAITLD  NORMAL  SI'RKADI  NO  I-TNCTION 


Assume  that  the  directional  spectrum  can  be  represented  by  the  form 


S(f,  )  =  S(f)  D( f ,  )  (5.1) 

where  S(f)  is  the  frequency  spectrum  of  the  sea  surface  and  D(f,  )  has  the  form 


D(f,  )  = 


.exp_[-L(  .1 


(5.2) 


where  and  are  functions  of  frequency.  This  functional  form  is  known  as  the 
wrapped  normal.  The  function  D(f,o)  has  the  following  Fourier  coefficients: 
exp(-n;>  y  / 2)  c o s ( n 0 <, ) 


(5.3) 


e x£( - rr_ f_  /2) 


sin(n  0) 


(5.4) 


It  can  be  shown  that  in  the  case  of  water  elevation  x  slope  and  y  slope  that 


C  ( f )  =  S(f)k2n  b?/2 


Sxx(f)  =  S(f)k?(l+^a,)/2 


Syy(f)  =  S(f )k? (1 -7  a? )/2 


(5.5) 

(5.6) 

(5.7) 


(Borgman , Hagan  ,  and  Kuik  1982) 

The  coherence  from  Equation  (4.17)  is 


C?  (f) 

v  w  '  ' 


Coh(f)  =  xyv 

s  TfT  s^TfT' 

'  yy '  / 


(5.8) 


xx 


Substituting  the  value  in  (5.5)  through  (5.7),  it  is  rewritten  to  be 


Coh(f)  =  —^20^.. 

exp(4o- )-cos  ( 2‘:  0 ) 

Thus  the  coherence  can  be  tabulated  for  varying  values  of  the  mean  direction 
and  the  spread  parameter  ■'  .  As  apparent  from  the  values  in  Table  7 
for  typical  values  of  o  and  ,  the  coherence  varies  from  0.0  to  around 


Having  this  range  of  values  for  Coh(f),  it  is  possible  to  develop 

tables  for  the  confidence  multipliers  tj  and  t?  from  Equations  (4.6) 

through  (4.15).  These  tables  appear  in  Appendix  C.  Appendix  D  contains  the 

corresponding  values  of  the  parameter  D  from  Equation  (4.3). 

The  confidence  interval  that  results  from  Equation  (4.7)  was  broken 

into  four  cases  depending  upon  the  values  of  t x  ,  t2  ,  and  S  ..(f) .  This 

xy 

makes  the  interpretation  of  these  intervals  difficult  since  cases  3  and  4 

result  in  intervals  of  infinite  length.  However,  another  look  at  Equation 

(4.7)  yields  the  confidence  interval  on  the  reciprocal  of  S  (f),  that  is, 

x  y 


(  ti 


,  t,>  ) 


(  svu(f)  s  (f)  ) 


for  S  (f)  0  and 

(  t;..  ■ 


(  .  t;  ,  ) 

(  swu(f)  sv  If)  ) 


for  S  (f)  0 

xy 


The  confidence  interval  on  1/S  (f)  is  more  easily  computed  and  exists 

in  a  meaningful  form  for  all  t  ;  ,  t  ,  and  S  (f)  except  S  (f)  =  0.0. 
Another  good  property  of  this  interval  is  the  fact  that  it  has  a  finite 


width  that  is 


_ J2 _ 

ti _ | 

S(f) 

3(f) 

xy 

xy'  '| 

1 

S  (f) 
xy  ' 

(tp  -  tj 

(5.12) 


since  t2  ■■  tj  for  all  tx  and  t2  satisfying  Equation  (4.7). 

The  relative  width  for  the  intervals  (5.10)  and  (5.11)  is  defined 

to  be 

RW  =  (t.  -  tj  (5.13) 

The  tables  in  Appendix  C  also  contain  the  corresponding  values  of  RW. 

In  Appendix  D  is  a  listing  of  the  computer  programs  that  compute  tj 
and  t?  along  with  a  description  of  how  they  can  be  made  operational  on 
any  computer  equipped  with  any  of  the  standard  FORTRAN  languages. 

The  confidence  intervals  on  1.0/Sxy(f)  for  the  values  in  Table  5 
appear  in  Table  8  along  with  W  ,  the  interval  width.  A  spread  parameter 
of  a  =  20  degrees  was  used  with  o  =  o  in  order  to  compute  coh(f)  and 
thus  tj  and  t2  •  Equations  (5.10),  (5.11),  and  (5.12)  were  used  to  obtain 
the  values  in  Table  8.  The  values  used  for  tj  and  t2  were  obtained  from 
Appendix  C. 


m 


Table  7 

Coherence  as  a  Function  of  the  Wrapped 
Normal  Parameters 


>  i  G I .  /i 


•;ta 

16.0 

16.0 

17.0 

18.0 

19.0 

2  0.0 

o . 

0.0000 

0.0000 

C  .0000 

c  .  0  0  0  0 

0.0000 

C  .  0  0  0  0 

C-  • 

.0152 

.0131 

.0114 

.0100 

.0087 

.0077 

4  . 

.  057  9 

.0603 

.04  39 

.  0  385 

.  03  3  9 

.0299 

6  . 

.1205 

.1056 

.0  92  9 

.  0620 

.0726 

.  06  44 

8  • 

.1541 

.17  19 

.16  2  5 

.1357 

.  120  9 

.1079 

1  0  . 

.  2705 

.  2422 

.2170 

.  1  946 

.17  47 

.  1570 

1  2  . 

•  3  440 

•  3113 

.  2  G  1  6 

.  2547 

.  2  3  04 

.  2085 

1  4  . 

.4113 

.  37  56 

.  34  30 

•  3129 

.  2  8  5  2 

.  2598 

16  . 

.4710 

.  4  34  1 

.  3995 

.  367  1 

.3370 

.  3  0  90 

1  6 . 

.  522  6 

.4655 

.450  1 

.4165 

.  3  84  7 

.  35  4  9 

2  0  . 

.  567  1 

.  5  302 

.  4947 

.  4605 

.  427  9 

•  3  96  8 

2  2  . 

.  6047 

.  5  6  G  6 

.  5  3  34 

.  4  992 

.  4662 

.  4  3  45 

2  4  . 

.  636  5 

.6014 

.566  8 

.  5329 

.  4  999 

.  467  9 

Table  8 

Confidence  Intervals  for 


k 

0 

coh 

20 

-9.1 

.13 

16 

-45.0 

.61 

18 

-15.6 

.29 

24 

6.9 

.08 

34 

12.0 

.21 

52 

5.9 

.06 

84 

4.5 

.03 

172 

-2.6 

.01 

Note:  Angles  in  degrees.  For 


^  -i  Confidence  interval 

ty  i.o/sxy(f> 


-.954 

(-3.34, 

1.18) 

-.299 

(-.66,  - 

-.03) 

-.115 

(-.32,  . 

05) 

.017 

(-.028, 

.066) 

.021 

(-.007, 

.052) 

.074 

(-.084, 

.240) 

.175 

(-.249, 

.611) 

-.315 

(-1.26, 

.62) 

;  =  20  degrees ,  1  -  »  =  .95. 


PART  6 :  SUMMARY 


A  method  for  computing  confidence  intervals  for  the  spectral  energy 
estimates  found  in  the  monthly  reports  of  the  Coastal  Data  Information 
Program  was  presented.  The  chi-square  values  presented  in  Table  1  enable 
the  quick  coni  station  of  confidence  intervals  for  any  location  and  time 
given  in  t h« *  >u 'ports. 

'he  p'S'P.go  1  i  tv  law  the  estimate  of  the  longshore  component  of 
rodiuti  -i  r.-  •  denoted  v  * '  was  found  to  depend  on  the  coherence  of 

•y 

too  lonosu  .'  '  .  e  •  -.  nt  with  the  Of f shore  slope  component.  Since 

•  ‘  bat.'.  ■  n>;  !•  ■  .  are  not  available,  either  finding  a  confidence 

1  nterva !  tn.it  g.  .  ft.  is  ai  the  distribution  of  the  spectral  density  estimates 

rather  tnar.  f  r .  •  !i  •  tr  i ;  u  f  i  an  .>?  S  (f)  or  obtaining  values  of  the  coherence 

•<y 

tv  assuming  a  parameter  i. ■:»><;  di  rectional  spectrum  must  be  considered.  The 
fir.*  approach  is  made  possible  by  assuming  a  value  of  the  apparent  angle 
Hid  then  nmi'ut.i  ng  tne  interval  of  Equation  (4.27).  The  second  approach 
re qui  res  the  assumpti on  of  a  wrapped  normal  spreading  function  (or  perhaps 
some  other  parameterized  directional  spectrum)  for  which  the  coherence  is 
derived  as  a  function  and  then  computed  for  specific  values  of  the  angle 
parameter  and  the  spread  parameter  (of  Equation  (5.22)  .  After 

the  coherence  is  computed,  the  values  of  t  and  t  are  obtained  using 
the  program  (i  1  Al .2  listed  in  Appendix  1).  Contained  in  Appendix  D  are 
listings  of  t  and  t  lor  various  values  of  the  degrees  of  freedom  k  , 
confidence  (oetfi;  i  ’nt  1-.  ,  and  inherence  COH,  The  first  method  ’provides 
intervals  on  the  quantity  s  ( f )  itself.  I  he  second  method  has  diffi- 
(iilties  with  cases  where1  the  constant,  f  is  less  than  zero  and  the 
constant  t  is  greater.  When  this  n.fiirs,  the  intervals  become  infinite 


in  width  and  are  therefore  difficult  to  interpret.  However,  the  confidence 
interval  for  the  reciprocal  of  S  (f)  does  not  exhibit  these  difficulties. 
It  is,  therefore,  presented  as  a  much  more  easily  understood  interval.  The 
interval  for  percent  energy  in  Equation  (2.12)  has  relative  width 

RW  -  k(l/Xkii/2  -  l/XkJ„i/2) 

[he  relative  width  has  values  varying  from  .36  for  period  band  6-4  to  1.4 
for  period  band  Id- 16.  RW  can  be  thought  of  as  the  resolution  of  the 
estimate  with  90  percent  certainty.  Thus  the  percent  energy  can  be  esti¬ 
mated  with  36  percent  to  140  percent  resolution  depending  upon  the  period 
banc  or  value  of  k  .  Accordingly,  the  interval  for  radiation  stress  has 
relative  width  RW  from  Equation  (5.13).  RW  is  listed  under  W  in  the 
table  of  Appendix  C.  For  typical  values  of  the  coherence  from  Table  8, 

RW  varies  from  2.11  to  5.90  resulting  in  a  resolution  of  211  percent  to 
590  percent  with  95  percent  certainty  on  the  estimate  of  1/S  (f)  .  In  all 

cases,  the  confidence  intervals  on  1/S  (f)  contained  the  value  zero; 

xy 

consequently,  inferences  about  S  (f)  are  difficult.  Due  to  the  apparent 

xy 

lack  of  resolution  and  to  the  difficulty  in  interpreting  the  intervals  on 
1/S  (f)  ,  it  is  concluded  that  the  method  given  by  Table  6  is  the  better 

of  the  two  tor  general  purposes. 
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APPENDIX 
PROOF  OF  CONDIT 


USING  RESIDUE 


Case  1:  Given  z  0  and  the  curve  is  chosen  to  be  t  =  Re  0 


In  order  to  apply  the  residue  theorem  it  must  be  true  that 


Qz( t)  M/'t 

for  t  1arge,  and  no  poles  of  Qz(t)  on  the  real  axis  (Marsden,  1973). 
Then  t ;  =  R  and 


Qz(t)  =  t-ia  (t-ib) 

has  poles  on  the  real  axis  only  if  a=0  or  b=0.  This  cannot  happen  due  to 
the  way  in  which  a  and  b  are  defined.  Now  consider 

t-  i  a  |  1 1-  i  b 1  =  I  Re  ^  - i  a :  |  Re  1  -  i  b !  ^  f  |  Re  ^ l*  |  -  |  i  a  |  |  1 j  Re  ^ '  \ 

=  ' R-(a; ; ’ R- I b ; |  =  j R: -R( ' a | + | b | )  +  |ab| 


then 


t-ia  t- ib  -  R  -R( ' a 1  + ; b ' )  +  | ab i  1 
and  if  R  is  large,  then  R  •  ja’  +  jb'  so 
t-ia  t-ib  R  -R(  a  +  b  ) ! 
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confidence  coefficient  =  1--.  from  Equation  (4.7);  D  =  parameter  from  Equation  (4. 
71  and  T2  =  t  and  t  from  Equation  (5.10);  and  W  =  relative  width  (t  -  t-). 
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=cent  city  array. 
i  IRS  1 


PERCENT  ENERGY  IN  3 AND 
(TOTAL  ENERGY  INCLUDES  RANGE  2043-4  5EC3) 


CAL 

3  I  G  HT 

TQT  EN 

BAND 

PERIOD 

LIMIT 

S  (SEC 

3 ) 

r  I  ME 

(  cm  ) 

(CM  30) 

22- 

22-18 

13- la 

16-14 

14-12 

12-10 

10-8 

8-6 

0220 

75  4 

-■=  <=  3 

1  5 

1.  7 

3  6 

8  0 

4  5 

20  0 

24  9 

21.  8 

0926 

91  1 

513  3 

1  5 

0  4 

3.  7 

9  4 

8  2 

20  5 

32.  5 

13,  6 

1  C  ">  1 

*  c.  . 

105  7 

693  1 

2.  2 

1  3 

1 .  6 

9  4 

18  1 

23  0 

20.  1 

10  9 

2 ' 20 

133  9 

1  120  3 

3  0 

0.  6 

0.  a 

3.  0 

43  8 

20  7 

12.  3 

7.  4 

0222 

118  6 

879  7 

3.  4 

0  5 

0  8 

2.  3 

20  0 

30.  3 

21  7 

10  1 

0R2C 

ICO  5 

630  7 

2  9 

0  4 

1  1 

3  3 

1C  1 

21.  2 

32.  0 

17  9 

1  52C 

102  9 

66  1  .  3 

1  9 

0  5 

1  5 

3  3 

10  7 

30  ? 

25  7 

14  5 

2119 

92  2 

531  2 

4  6 

0  3 

1  6 

2.  3 

5  9 

23.  4 

27  1 

14  4 

0221 

88.  5 

46q  3 

4  3 

0  3 

2.  5 

1.  8 

2.  R 

21  3 

25.  6 

17  9 

0920 

75.  5 

356.  2 

4  3 

1.  3 

1 .  6 

2.  4 

—  — 

7  7 

42.  6 

23.  0 

1519 

75.  3 

354.  6 

7  6 

2.  0 

5  2 

2.  R 

1 .  9 

11  2 

35  7 

19  4 

2120 

62.  4 

243.  0 

7  2 

1  5 

16.  1 

3.  6 

4  6 

11.3 

26  0 

15  9 

0220 

81.6 

416.  2 

9  3 

0.  8 

3.  0 

2.  7 

2  3 

4  0 

43.  7 

24  6 

0R22 

72.  5 

323.  4 

9  1 

1.  0 

4  3 

4.  4 

3.  0 

7.  9 

31.  6 

21  5 

1  524 

75.  1 

352.  0 

6.  3 

1.  3 

4  7 

4.  2 

3.  9 

6.  7 

42.  8 

16.  4 

*11  •’*! 
c_  *  e. 

70.  6 

311.9 

8.  7 

6.  3 

5.  7 

6.  7 

4  1 

9  9 

25  3 

10.  9 

0321 

718 

322.  1 

12.  0 

2.  3 

2.  7 

2.  9 

6  4 

20  4 

29  4 

13.  6 

1C07 

69  0 

297  5 

11.0 

3.  8 

3.  7 

8.  5 

2.  3 

17.  1 

25.  3 

19.  3 

1  523 

61.  7 

233.  0 

9  3 

3.  3 

3.  5 

6.  8 

3.  0 

14  5 

30  5 

15  5 

2120 

60.  1 

225  5 

8.  a 

3.  1 

8.  1 

6.  6 

4  a 

13.  3 

27  2 

14.  0 

0320 

59  7 

222.  9 

14.  6 

4  5 

5.  4 

7.  6 

2.  7 

7  4 

35.  4 

15.  5 

0925 

61.  6 

237.  0 

6  9 

20.  0 

8.  3 

5.  3 

2.  3 

4  0 

35.  7 

9.  0 

1  524 

57.  3 

205.  4 

8.  1 

22.  1 

10  2 

9.  2 

3.  8 

5.  1 

20.  9 

12.  4 

2121 

58  0 

210  0 

5  3 

30  4 

10  5 

6.  2 

3.  7 

3.  6 

16.  4 

10  7 

0320 

54  0 

182.  5 

13  1 

24  0 

7  9 

8  6 

7.  5 

4  6 

110 

12.  9 

0923 

64  0 

255.  9 

4  2 

21.  7 

29  2 

4.  7 

4  0 

2.  3 

17  6 

7.  5 

1522 

47  5 

141.  2 

2.  2 

10  6 

37  1 

16.  8 

4  4 

4  1 

10  0 

8,  4 

2121 

53.  1 

176.  1 

4.  2 

13.  8 

37  3 

10.  6 

7  3 

7.  3 

116 

4  7 

0320 

51.  7 

167.  2 

10.  0 

10  0 

24  7 

18.  7 

6.  2 

7.  3 

11.5 

6.  1 

0R22 

74  2 

344  4 

1.  3 

5  1 

49  3 

17.  5 

2.  6 

1.  1 

16  9 

3.  6 

1  520 

51  9 

168  6 

4.  3 

6.  3 

24  a 

35  R 

9  6 

5.  7 

5.  5 

5  9 

2122 

81.  3 

412.  8 

1.  4 

9  2 

15.  7 

14  6 

3.  6 

2.  1 

10  6 

25  5 

0320 

78.  0 

300.  3 

2.  2 

3.  2 

10  1 

27.  6 

13.  6 

2.  4 

17  3 

15.  3 
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APPENDIX  B 

PAGES  FROM  ORIGINAL  REPORT  FOR  WHICH  EXAMPLE 
CONFIDENCE  INTERVALS  ARE  COMPUTED 


and  if  R  is  large  enough,  then  R-(  a  +  b  )j  '•  1.0  and  thus  for  sufficiently 
1 arge  j  t !  =  R 

t-i a | ! t-i b  •  R 

There  fore 

|  Qz  ( t ) !  =  f  [  t-i  a  |  j  t-i  b ;  >  R"  ' 

and  then 

lQz(t) !  <  1/R 

but 

1/R  y  1/R  =  1  / 1 1 f 

The  boundedness  condition  is  therefore  satisfied  by  M  =  1.0. 

Case  2:  Given  z  0  and  the  curve  F  is  chosen  to  be  t  =  Rel!  0  £  '■  <  -. 
Then 

!  t- i a !  |  t-i b |  =  [  Re ^  - i  a  |  |  Re 1 '  ~ i b j  >  |  R- |  a  |  |  j  R- j  b|  | 
and  the  result  will  be  identical  to  that  of  Case  1. 
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APPENDIX  D 

TABLE  GENERATING  PROGRAM 


iV.cn  (it  1  on  of  Proyrai"  CFAB2 


r  J\  C  ij  t< » .  V  I  i\  r  c!  (  I  .  ;  '  V  ,  J  U  VP  u  1  ,  T  A  P  V  V  ) 

D  j  i;F!. ..  i  0 N  T  A :  (  9 1  j  ,  ■_  ),:;!.'(  5  ) 

DATA  (  !<:;  (  I  )  ,  1=  1  ,  ,}  /  AD  .  ,  1  0  .  ,  8  .  ,  9  -  : 

A  -  •  C  V 

CC -  1  .  -  A 

!  I T  =  9  c 

rnt  =  :;t 

D  T  -  .  9  /  i-  :i  T 
DC  DO  K  =  1  ,  9 
PRINT*  ,  RIK  K  ) 

H  K  =  2  *  h  .‘i  (  K  ) 

DC  10  1=1  ,  N  T 
C  C I !  = ( I  )*DT 

DEL  =  Mw  AT  (  CCH/  (  1  .  -CCli )  ) 

D  =  1  .  /  D  e;  1. 

C  A  L  l,  Z  i.  ARC  K  (  L  E  L  ,  iu.'  (  K  )  ,  A  ,  C  L  ,  C  U  ) 

TAD  (  I  ,  1  )  =COii 

TAD ( I ,2 ) =D 

TALK  I  ,  j  )  =CL*D/Rt:  (  K  ) 

talk  if4)  =cu  * l/rn  (  k  ) 

TAP  (1,5)  =  T  A  Li  (  1,4)  -TALK  I  ,  3  ) 
CONTINUE 
W  R  I  T  E  (  2 , 9  1  ) 

WRITE ( 2 , 92  )  RK 
WRITE ( 2,95)  CC 
W  R  I T  E  (  2  ,  9  3  ) 

DO  30  1=1,30 

W  R I T  E  (  2 , 9  4  )  (  T  A  P  (  I  ,  J  )  ,  J  =  1  ,  5  ) 

W  R  J  T  E  (  2  ,  9  1  ) 

WRITE ( 2 , 92  )  RK 
WRITE ( 2  ,  95 ) CC 
UR I T E ( 2 , 9  3 ) 

DO  40  1=31,60 

WRITE(2,94)(TA1KI,J),J=1,5) 

W  R  1 T  E  (  2 , 9  1  ) 

WRITE ( 2 , 92 ) RK 
WRITE (2,  95;  CC 
W  R I  T  E  (  2 , 9  3  ) 

DO  50  1=61,90 

W  n  I  T  E  (  2 , 9  4  )  (  T  A  IK  I  ,  J  )  ,  J  =  1  ,  5  ) 
PRINT*,  " R  K  =  " , R  K 
CON  TI NU  E 
FCRi.AT  (  "  "  ) 

F  0  r;  II A  T  (  7  X  ,  "  C  0  N  F  1 1;  E  N  C  E  MULTIPLIERS 
£  F  5 . 1  , "  DEGREES  OF  FREEDOM" ) 

;  F  0  R 1 1 A  T  (  7  X  ,  "C0li",7X,  "D",  <?X,  "TI  ",  SX 
FORMAT!  IX, 5 F  10.4 ) 

F  C  R  M  A  T  (  7  X  ,  "  C  0  N  F  I  D  F.  II  C  F.  C  0  E  F  F  I  C  1 1. T 


•CCli  )  ) 


»  !•'  n:.-:  ll  a p ci; (  :>  ,  r l , ,  cl  ,  it  ) 

■  1  •  ;  !  t  ■  c  i;  c  c  c  cccccctccccccccc  c  c  c  c  o  c  c  c  c  c  c  c  c  c  c c  c  c  r;  c  c  c  c  c  c 

'  tei.;.  .t::;  .ujtink  gkaru-ec  for  :•  tap  ting  FuIntl  fc;:  c 

TL!  END  I. C- 1:  T I !.  K  13  i  CF.CT .  C 

D=  NON  C  h li  T i;  A  L  I  T  Y  i* A  R  A  I!  1.1  L  K  C 

1. U  =  11  EC  REE  3  FREED  Cl!  C 

c  a=3iun  ;fic  a::ch  level  c 

CL.  =  L C •  ER  CON  F  L  DLL  Cil  L  Til  T  (J 

r  C  L  =  U  r  P  E  F.  C  0  F  3  DELI  Cel  LLIUT  C 

C  r 

c  c.kai-C:'  compute:;  a/.  ;.:id  i-a/?  t hen  r  i  r; d c;  pc  in  to  c 

c  THAT  li  A  V  L  VALLE.-.  OF  THE  FUNCTION  F  E-TAR  THAT  C 

C  E  T  A  I  I:  L  I.  A/C  AN  l'  1  -A/C  .  BISECT  13  THEN  CALLED  C 

C  TWICE  TC  HERO  IN  ON  THE  LOIHIE  ALL  UPPER  CON  II  DEUCE  C 

c  L  It:  I  TC .  !.’ Y  11.  ANL'EK’.:  UNIV.  OF  V.  Y .  C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

I T  C  =  0 
A  2  -i\  /  2  .  0 
P-l  .-A? 

IF(  Rt;  .  Lh.  0  .  )  PRINT6,  "ILLEGAL  i:L  " 

I  F  (  N  .  L  L .  0  .  )  R  E  T  0  R  N 
T  =  I)  *  R  :: 

10  F=F3TAR( T  ,  D, RL ) 

I  F  (  F  .  L  E .  A 2  .  OR  .  F  .  G  E  .  A  3  )  G  0  TO  L  0 
T  U  =  T 
TL  =  T 

R  1  N  C  =  (  .  1  )  »  30  RT  (  (  D  *  6  A  +  J  .  C  )  *  FI  1;  ) 

2  0  T  L  3  =  T  L 

T  L  -  T  L  -  R  3  IJ  C 
F  L  =  F  3  T  A  V  (  I  L  ,  D  ,  li  N  ) 

I  F(  FL.  L T. A 2 ) GC  TO  3  0 
CO  TO  0 
J  0  T  U  3  =  T  U 

T  U  =  T  U  +  f;  1  N  c 
FU  =  FETAL  (TU,D,  RM  ) 
i  F( FU . GT. A3 ) CO  TO  A 0 
GO  TO  30 

*i(  CALL  Bf  3ECT  (  TL  ,  TL3  ,  CL  ,  AP  ,  D  ,  RN  ) 

CAL  I.  f.I  3LCT(  T 03, 7U  ,  CU  ,  A  ,  D  ,  RN  ) 


(j  I  F  (  F  .  L  L.  A,,.  )  T  =  T+K i  N  C 
it-  (  F.  OF.  A  •'  )  T  =  T-niNC 
i  T  C  =  i  T  C  +  1 

i  r  C  !  TC  .  GT.  1  )  PR  1  NT*,  "TL  Oil  LI."  ,  17C 


s u  l  i-o  utile  ;  i  s E  c t  (  t  l  ,  v u ,  v ,  f  (  >  s >  R ’  ) 

C  c  eccccccccccccccccccccccccceccccccccccccccccccccccc c c c c c 


b j  j e c t  uses  the  ^section  r:i; th of-  to  zero  in  0::  c 

A  VALUE  T  THAT  CORRESPONDS  WiTI!  A  VALUE  CF  TiiE  C 

FUIICTIOL  F  S  T  A  R  (  T  )  =  T  0  .  C 

V  A  F  1  A  L*  I.  ES  :  C 

TL=LOWER  GUESS  Af  T  PASSED  FROM  C 

SUE  ROUTINE  SEARCH  C 

TUsUPPER  GUESS  AT  T  PASSED  FLO!.  C 

SULROUTINE  SEARCF  G 

T  =  F I N A  L  APPROXIMATION  TO  T  G 

F0= VALUE  OF  FUNCTION  FSTAP  THAT  G 

FSTAR(T)  MUST  EE  CLOSE  TO  C 

D=  EG!;  CFETR  ALITY  PARAMETER  C 

RN=DEGREES  FREEDOM  G 

is  Y  !1.  ANDREW  UtllV.  OF  UY.  G 


vccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

TOL  =  .  000  0  1 

ST  D  =  SQ  i'T  (  (  D  *■  *  2  +2 . 0  )  *  R li ) 

10  I F ( ( AES( TU-TL) /STD ) . LT. TOL)  GO  TO  2 0 
TLS=TL 

T  L  = ( TL+TU) /2.0 
F  =  F S T A R  (  T L  ,  D,  RN  ) 

IF( F. LT. FO  )  GO  TO  10 
TU  =  TL 
TL=TLS 
GO  TO  10 

20  T=(TL+TU)/2.0 
R  E  T  U  R  f  J 

EL  D 


r ) 

l  C  i :  C  ci.  c  C  (.:  C  C  i.  L  c  0  G  G€  C  C 

!  giFruieg  value  .•  ;  (  v;,: 

:>  i  0  A!  i  V  A  i  Ini-I.r.  ;  TI;  C  i  A  A  A  C  T  L  ft  I  G  T  J  C 

Vi  or-  o (  t  )  =  i  / (  i  -l  *::«T-rT*£i  )  r  •  nr:;. 

!  T  (  -  1  .  0  )  . 

A  V  L  i:  G  : 


=  cpp.  i:  LJi’jT  or  integration  c 

b-  tier;  gen  tr  a  l.  it  y  par  a i : i; t l: i-:  c 

HI.  U  =  DEG  HKES  frf.ecc::  c 

uses  the  fact  that  integrals  or  the  g a : : ;i a  c 

CAN  LK  WRITTEN  AG  PCISSOIi  SUES  (  GEL  Cii.IV).  C 

.  ANDREA  UN  IV.  OF  Y  .  C 


CCCCL’CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC'L 

.  0 

I; u 

+  F  i  (  L  *■'  A  +  -4  .  )  )  /  2  •  0 
L  -  GO  I T  (  D  8  A  +  L  .))  /2.0 
-■  FT  (  L  -  *  E  +  A  .  ) 

R  =  0  .  0 

.11.0.0)00  TO  ij  0 

A  1  =  1,  Nl!  I. 

1 


u  +  x  j 
u  -  v.  j 

* 

i-  /  ( ii  i  ' *  i  r  ;f)  +  (  i .  o-poxGS(  idif)  )  /  (  a:  g(  :  r  } * * ;  r ;  r ) 

=  ;  -  :  A  ;•’  +  1:1.  *1:  I  /  (  R  j  o  *  I  G u ;  i ) 

J(  +  )/(  xj+i . ) 


END 


FILMED 
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